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ABSTRACT 

We consider the formation and migration of protoplanetary embryos in disks around the stars 
in tight binary systems (separations ~ 20 At/). In such systems, the initial stages of runaway 
embryo formation are expected to only take pl ace within some cr itical disk radius Ucru, due to 
the perturbing effect of the binary companions (iThebault et al.l2009l) . We perform n-body sim- 
ulations of the evolution of such a population of inner-disk embryos surrounded by an outer- 
disk of smaller planetesimals. Taking Alpha Centauri-B as our fiducial reference example in 
which flcri, ~ 0.7 At/, and using a Minimum Mass Nebular Model with 2 oc a'^l^, we find 
that within IQ^yrs {\(fyrs), systems will on average contain embryos which have migrated 
out to 0.9 AU (1.2 At/), with the average outer-most body having a mass of 0.2Me (0.4M®). 
Changes to increase the surface density of solids or to use a flatter profile both produce in- 
creased embryo migration and growth. At a given time, the relative change in semi-major axis 
of the outer-most embryo in these simulations is found to be essentially independent of Ocnv, 
and we note that little further embryo migration takes place beyond 10^ yrs. We conclude that 
the suppression of runaway growth outside Qcrit does not mean that the habitable zones in 
such tight binary systems will be devoid of detectable, terrestrial mass planets, even if acru 
lies significantly interior to the inner edge of the habitable zone. 



1 INTRODUCTION 



in multiple systems 
20% of detected exo- 



The majority of stars occur 
('Eggleto n & TokoviniiJ 12008*) and ~ 
plane ts have been found in binary systems (Desidera & Barbieri, 
l2007h . hence the importance of studying planet formation in such 
systems. Most of these systems have very large separations for 
which the influence of the companion star on the planet region 
is probably limited. However, 3 planets have been detected in 
tighter binaries whose separation, at, is around 2Q AU, with 
one such system (y Cephei) having a planet wi t h a m inimum 
mass of 1.6Mj„,, orbiting at 2.1 At/ iHatzes et alj j2003h . While 
such systems might harbour planetary orbits which have a long 
term stability in regions up to 3 to 5 A (7 from the primary (E.g. 
iHolman & Wieger3ll999l) . they are thought to present much more 
chall enging environments for planet formation in these regions 
(E.g. lKlev & Nelsoi] ( 12007 ^ and references therein). 

In recent years, this issue of planet formation in tight bina- 
riesQhas been investigated in numerous studies. They have shown 
that the late stages of planetary accretion, i.e. the final assembly of 
large Lunar-sized embryos, can proceed in a domain almost as large 
as the orbital stability region ( Barbi eri et alj|2003 : lOuintana et al.l 
I2OO7I : [Ouedes et al. 2008). However, the earlier stages, starting 
from kilometre-sized planetesimals and leading through runaway 



' We refer here to "tight" binaries as systems with separations around 
20 AU. Systems with even smaller separations are left out of this study since 
they usually imply unstable orbits in most of the potential terrestrial planet 
region and are thus of limited interest for the present problem 



and oligarchic growth to the embryos, are much more affected by 
the binary companion. Indeed, the coupled effects of secular per- 
turbations from the companion star and gas friction lead to strong 
orbital phasing as a function of planetesimal sizes. Such phasing 
results in high impact velocities between differently sized objects, 
causing accretion to be significantly slowed or even turned to frag- 
mentation for a large range of parameter space. As a result, only 
the innermost regions of the protoplanetary disk, beyond a critical 
distance acrit from the central star, might allow the in-situ forma- 
tion of planets (E^. lThebault et alj200a : |Paardekooper et alj200^: 
Xie & Zhou 200g). For binaries with a^ ~ 20 AU, a„,f could be 
(depending on the system's ecce ntricity and mass ratio ) well inside 
1 AU (see for example Fig.8 of IThebault et ai]|2006l) . This could 
pose serious problems for explaining the presence of a planet like 
that of the y Cephei s ystem at 2.1 AU (see for instance the numer- 
ical investigations of .Paardekooper et al. 2008; Xie & Zhou 2008b . 
or that at 2.6 AU in HD19688 5 (although the orbi t of the binary is 
here not fully constrained, see lCorreia et al.|[2008l) . 

A good illustration of these problematic issues is the case of 
our closest stellar neighbour, the tight binary Alpha Centauri {a 
Cen). To date, no planets have been found in the system jEndl et al] 
( I2OOII) r ule out any planets w ith masses > 2.5Mj„,,), but a recent 
study bv lOuedes et alj ( I2OO8I) suggests that, ;/ terrestrial planets do 
exist in this system, then a Cen-B would be the ideal candidate for 
the detection of potential terrestrial planets via the radial velocity 
technique, due to the proximity and quiet nature of the star. How- 
ever, despite long-term stability as well as embryo-assembly stud- 
ies which have given very promising results for the regions interior 
to 2Af/ dHolman et alJI 19971 ; Isarbieri et al.ir2002l : lOuintana et al.l 
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l2002l : lGuedes et al.l2008t) , lThebault et"ai] j2008ll2009l) have shown 
that the 0^™ for planetesimal accretion could be located as close 
as 0.5 AU from each star. This would place it at the inner edge 
of the habitab le zone (HZ) of a Cen-B (~ 0.5 - 0.9 AU, see 
iGuedes et all20 08). but well inside that of a Cen-A (~ 1.0-1.3 AU 
seelBarbieri et al.,.2002) . Although a recent study by IXie & Zhoij 
shows that a small inclination between the planetesimal disc 
and the binary orbital plane would help push the critical radius 
further out, the in-situ forma tion of planets i n both HZs remains 
highly problematic. However. iThebault et alj lioOSl) point out that 
the presence of planets in regions beyond a„ir cannot be ruled out if 
additional mechanisms are invoked. Such mechanisms include (i) 
a change in dynamics after gas dispersal, (ii) a change in the bi- 
nary orbit properties at some early point in the history of a Cen, 
or (iii) outward migration o f the planetary embry os which formed 
in regions with a < 0.5 At/. IThebault et all ( 12009 ^ go on to discuss 
points (i) and (ii) in some detail, ruling out the former and acknowl- 
edging the latter as a possible solution, but leave the detailed study 
of (iii) to future investigations. 

Our aim in this paper is to address this crucial issue in more 
details. The main question we want to answer is the following: 
could large embryos, accreted in the inner regions interior to a„ii, 
later migrate outward and end up beyond flcnr? And if so, could this 
radial displacement be enough to place planets in the HZs of tight 
binary systems? 

We note that the migration of planetary embryos has been 
observed in ma ny previou s inves tigations. Among the many ex- 
amples of this, IChamberj ( 1200 ll) demonstrate mixing and rear- 
rangement of emb ryos due to embryo-embryo scattering, whilst 
iKirsh et all ( |2009|) show embryo migration which is driven by 
embryo-planetesimal interactions. We would thus very generally 
expect to observe embryo migration within our simulations and 
simply wish to understand the degree and nature of the phe- 
nomenon. 

We numerically investigate this problem, using an N-body 
code, taking the a Cen system as a typical reference example of 
a tight binary. 

2 MODEL AND METHOD 
2.1 The Disk 

We start by considering the protoplanetary disk around one star in a 
binary system, just after the dissipation of the disk gas. The disk is 
thus taken to consist solely of solid material distributed with a sur- 
face density profile S = Soa"", where is the normalised surface 
density at lAU (in g cm"2), a is the semi-major axis (in AU) and 
Q- is a positive constant. Unless otherwise stated, we take Zq = 10 
and a = 3 /2 as used in the nominal Minimum Mass Nebu lar Model 
(MMNM. Eavashil ( Il98ll) ) case of .Thebault et alj ( l2009h . 

We divide the disk into two zones, (i) an inner zone, a < a„i,, 
where runaway growth is assumed to have taken place and embryos 
have grown to isolation, and (ii) an outer zone, a > a„ii, where run- 
away growth was suppressed and the material in the disk is com- 
posed of smaller mass planetesimals. We take a Cen-B as our refer- 
ence case in the majority of our simulations, taking a„,r = OJ AU 
for such a system. We then set inner and outer simulation bound- 
aries at 0.2 and 1.5 AU respectively, meaning that the total mass of 
embryos in the inner zone is ~ 2Me, while the total mass of the 
outer zone is also ~ 2Me (the width of the outer region was chosen 
to give a total mass approximately equal to that of the inner zone 
when cicrit = OJ AU). 



2.2 The Inner Zone 

We expect that in the regions of the disk interior to 0.7 AU, where 
runaway growth can take place, growth to isolation will have oc- 
curred prior to the gas disk dissipation, resulting in a chain of iso- 
lated embryos. 

To check the validity of this assumption, consider the growth 
timescale of a protoplanetary embryo of mass M embedded in a 
population of pla netesimals of mass m in the absence of an external 
perturber: this is jKokubo & Id3l2002h : 

UOgcm-2; \m^) llAU/ \mJ 

^i 2400gcm-^ J lloT^) 

where E(3n,o is the surface density of gas in the disk around the 
star of mass and K is an integer used to specify the width of 
the annular feeding zone, Aa = Kyh, in terms of the Hill radius, 
yh = a(M/3My,)'''. Such an embryo would thus grow to an isola- 
tion mass 

M,,„ = (2nKI.o)"HiM,)-"'-\^-^j (AU? (2) 
in a time, t ~ 3t, where the factor of 3 arises on integrating t oc 

This calculation ignores the refinements of the IThebault et al.l 

( l2009t) work on the suppression of runaway growth in the presence 
of a perturber. Critically, Eqn[T]requires that embryos and planetes- 
imals of mass ~ 10'**g (~ 5km) have already formed and continue 
to exist. Despite these obvious drawbacks, we employ Eqn[T]to give 
us an approximate estimate of the timescale for growth to isolation, 
and suggest that because (a) t is less than lO^yrs at small semi- 
major axes, (b) a putative gas disk around Alpha Cen is likely to 
have dispersed on a > 10' year timescale and (c) the suppression 
of the growth rate in the inner region is unlikely to be more than 
an order of magnitude, it thus seems likely that growth to isolation 
will have occurred in the region interior to 0.7 AU before the gas 
disk dissipates. Hence we feel it is reasonable to assume that there 
exists a chain of isolated embryos at a < a„„. It is unlikely that 
any significant evolution beyond this stage would have taken place 
at the time of disk dispersion, due to the very long orbit-crossing 
timescale in the presence of gas damping. 

For K = 7, Ho = Wgcm^^ and a ~ 0.5 AU, the isola- 
tion mass will be Mj^o ~ 2 x 10^''g and the isolation separation 
Aa,j„ ~ 0.02 A (/, meaning that e < 0.02 for non-overlapping or- 
bits. As such, a reasonable set of initial conditions for the inner 
regions of our simulations seems likely to be a chain of isolation 
mass embryos, around 50 in number, with average eccentricities 
<e>~ 0.01. 

We simulate these isolation mass embryos as a population 
of fully interacting massive particles within the n-body code, i.e. 
the embryos interact gravitationally both with themselves (embryo- 
embryo) and with the smaller planetesimals (embryo-planetesimal) 
which occupy the outer zone. Note that in the majority of the sim- 
ulations there are no planetesimals within the inner zone. 



2.3 The Outer Zone 

The mass of solid material in the outer disk (~ 2M^) is divided 
into a population of equal mass planetesimals, which are randomly 
distributed radially according to the a"^'^ profile, and the majority 
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of the simulations use 100 planetesimals in the outer zone, meaning 
that each has a mass m ~ 0.02Mg,. 

The planetesimals are modelled within the n-body code as test 
particles with an associated mass (massive test particles), meaning 
that they interact gravitationally with the embryos in the system, but 
they do not undergo planetesimal-planetesimal scattering amongst 
themselves. Some details of the effect of this approximation are 
given in 93.51 



2.4 Numerical Method 

We use the mercury n-body package of Ichambersi ( Il999h to per- 
form simulations of the evolution of such a protoplanetary disk with 
the aim of understanding the degree to which embryos can migrate 
from the inner zone (a < a„it) to the outer zone (a > an-n). 

A summary of all the simulations performed is given in Table 
[T] Unless otherwise detailed in that table, the simulations were set 
up and performed as follows: 

• As the gas disk is assumed to have dissipated prior to the start 
of our simulations we do not impose gas drag / damping on any of 
the simulated embryos or planetesimals. 

• The solid disk is split into two regions: (i) The inner region, 
0.2 < a < acrif, populated with isolation-mass embryos, and (ii) 
The outer region, a„j, < a < 1.5, populated with a swarm of lower- 
mass planetesimals. We take a^n = 0.7 AU in the majority of our 
simulations, 

• The majority of the simulations are based around a single star 
and do not include the binary perturber in the n-body simulation. 
We note that, perhaps surprisingly, excluding the binary perturber 
from the simulations does not strongly affect the evolution in the 
semi-major axis, mass or eccentricity of the outward-migrating em- 
bryos and hence we omit the binary perturber from the majority of 
our simulations. We delay a more detailed examination and discus- 
sion of this result until 93.5.21 

• We run mercury in "hybrid" mode for almost all simulations, 
allowing rapid simulation via symplectic integrator when there are 
no close encounters, but also allowing accurate resolution of any 
close encounters (by using the Bulirsch-Stoer integrator). However, 
for any simulations which include the binary companion, we solely 
use the Bulirsch-Stoer integrator. 

• Collisions occur when two bodies approach to within the sum 
of their radii, with the radii being calculated assuming a density, p = 
3 g cm"^^. Any collisions are assumed to result in perfect mergers. 

• For any given parameter set we conduct 10 version of each 
simulation, keeping the number of embryos and planetesimals con- 
stant across all 10 simulations, but randomising the orbital angu- 
lar elements of these bodies. The results of such an ensemble of 
simulations allows some insight into the variations inherent in the 
chaotic evolution of the system. 

• In our fiducial simulation (set A), we have ~ 56 embryos in the 
inner zone and 100 planetesimals in the outer zone. Generally, we 
have kept these numbers constant in the other simulation sets. How- 
ever, variations in other parameters (E.g. surface density) can indi- 
rectly result in the initial number of embryos varying (as a result 
of Eqn.|2}: where this occurs we note the change and also conduct 
additional simulations to check that the change in embryo number 
is not affecting the results. More details on the numbers of embryos 
and planetesimals in each simulation set are included in the perti- 
nent sections. 



3 RESULTS 

We provide a brief summary of the simulations performed and the 
key results obtained in Table [T] 



3.1 Typical Results 

Consider a MMNM disk in which isolation mass embryos are ini- 
tially confined to the inner region (0.2 - OJAU), while the outer 
region (0.7 - 1.5 AU) has had runaway growth suppressed and is 
populated by smaller mass planetesimals (simulated using massive 
test particles). We conduct a set of 10 simulations of such a disk, 
assigning to each simulation a different set of randomised angular 
orbital elements. We label this fiducial set. A, and display in Fig[T] 
the full orbital evolution of all of the particles in one of the simula- 
tions. 

We find that at lO^yrs, the initial 56 embryos have reduced 
(via collision) to ~ 10 in number, with the outer body having mi- 
grated outwards to ~ 0.9 AU. If we define the relative change in 
semi-major axis, 

^omer.final ^omerJniTm! 

r} = , (3) 

^ouTerJnitial 

then we find that rj = 0.29!° jJ^, where the errors indicate the varia- 
tions to the upper and lower quartiles. The mass of this outer-most 
body has grown from 0.06 to 0.2Me, primarily through the accre- 
tion of planetesimals. 

Continuing our simulations of set A onto lO'yrs (set B and 
Fig[T), we find that on average: (i) Further collisions have taken 
place, reducing the number of embryos within 1.5Ai7 to ~ 6 and 
(ii) The outer-most embryo has now migrated to ~ 1.25 A f/ (giving 
77 ~ 0.8^5524), while its mass has further doubled to ~ 0.4MiB. 

Only a very small amount of additional migration is observed 
between 10^ and 10^ yrs, suggesting that the degree of migration 
observed by lO' yrs suffices to give a reasonable approximation 
to the overall degree of migration expected. We should emphasise 
though that by this stage of the simulations the outer-most em- 
bryo has migrated to a position rather close to the outer edge of 
the simulated planetesimal zone at 1.5 At/, so any future simula- 
tions which wanted to properly examine migration to 10** yrs and 
beyond, would have to be run with many more planetesimal bodies 
spread out to much larger semi-major axes. 

We thus find that the embryos which start with semi-major 
axes a < a„ii will migrate outwards over time as they interact with 
the exterior planetesimals. This means that despite the early sup- 
pression of runaway growth in regions beyond a„,7, the Habitable 
Zone (HZ) of the disk (which would initially have been devoid of 
large bodies if «„,, lay at a smaller semi-major axis than the inner 
edge of the HZ) can become populated with large protoplanetary 
embryos within the relatively short time of 10^ years. 

We note that the masses of these outer bodies in our simula- 
tions remain rather small compared to t he detection limits of cur- 
rent technology such as Kepler (see E.g.'Deegl'2002')), but we note 
that at the 10' - 10** year point in our simulations, the systems 
are still rather densely packed (6 bodies within 1.5 At/). We ex- 
pect that extended simulations (to 10' years) would suffer further 
coUisional evolution, yielding systems with rather fewer but more 
massive bodies, bringing many within detectable limits. 
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Table 1. List of Simulations Performed. Note that refers to the surface density of embryos, whilst Z/j refers to the surface density of planetesimals. The 
migration of the outermost embryo is expressed via the quantity = (Ooufer-/mn/ - <^ouier,mitial)/Oouter,imiial- 



Set 



Parameter Variation 



Domain of Variation 



Key Results 



Set A Fiducial MMNM simulation 

Set B Length of time over which Simula- 10^ — 10** yrs 

tions are conducted, T 



Set C Surface density normalisation of 5 - 50g cm"- 

(fully interacting) embryos in the 
inner zone, l-Ejimer 

and ratio of outer surface den- 0-3 
sity to inner surface density, 

^PJ nner I^EJnner 

Set D Surface density profile of the disk, -1.5 —> 



SetE Position of the Inner-Outer Bound- 0.4- 1.0 AC/ 

ary, a„it 

Set F Number of (fully interacting) em- — 100 

bryos in the inner zone, NE.inner 
and Eccentricity of (fully interact- 0.01 — 0.5 

ing) embryos in the inner zone 

^EJnner 

Set G Fraction of the inner disk which is 0. 1 — 1 

composed of embryos, / 



SetH Number of (massive test-particle) 100- 1,000 

planetesimals in the outer zone. 

Set 1 Fully interacting planetesimals in 100 - 250 

the outer zone, NpEuU.Omer 
Set J Addition of binary perturber to the 

simulated systems 



Outer most embryo migrates to 0.9 AC/ (;; = 0.3) within 10*jt.s and has 
a mass of 0.2Me . 

As we lengthen the simulations by a factor of 10 the relative change in 
semi-major axis more than doubles (;;(; = 10*) = 0.3, ri{t = lO') = 
0.8) and masses almost double. Little outward migration is observed 
between lO' and 10** years. 

Increased surface densities increase outward migration and cause order- 
of-magnitude increases in masses of the outer bodies. 

Higher ratios of planetesimal to embryo surface density cause increased 
migration and order-of-magnitude increases in masses of the outer bod- 
ies 

Shallower surface density profiles cause increased migration due to the 
relatively larger amount of exterior mass, with = 1.0 at yrs for 
0- = 0. 

There is a large scatter in the degree of relative migration, but are ap- 
proximately consistent with being independent of cir„7 ■ 
Both the amount of outward scatter and the final mass of the migratory 
planet are strongly dependent on NE.inner- 

The initial eccentricity is unimportant for the migration as long as 
^E.iimer < 0.5, but does play a role in determining final masses. 

If the total suiface density is held constant, then lower / results in 
slightly more migration, but slightly less massive planets. 
If the total surface density increases at constant density of embryos, then 
both outward migration and mass are increased. 

When the individual mass of the planetesimals is reduced, an insignifi- 
cant decrease in the amount of migration is observed. 

The use of massive test particles as opposed to fully-interacting particles 
for the planetesimals does not strongly change the degree of migration. 
The addition of a binary perturber leaves migration distance almost 
identical and slightly decreases mass growth. 



3.2 Variation in Surface Density 

We now look at the results for set C in which we vary the surface 
density of the inner zone between 5 and 50g cm"^, and also look 
at the effect of varying the surface density ratio between the inner 
and outer zones (i.e. a scenario in which the inner zone is depleted 
compared to the outer region, or a scenario in which the outer re- 
gion has become enriched). Note that the simulations with surface 
density ratios of 2 or 3 were extended outwards so that the outer 
edge of the simulation zone was at 3.5AU, rather than 1.5 AU, to 
allow for the possibility of significantly greater outward migration 
(the number and total mass of the planetesimals was increased ac- 
cording to the surface density profile). 

We plot in Fig |2] the results of the surface density variations 
after 10* years. We see that in general, moving to a higher absolute 
surface density increases the degree of outward migration, as does 
a move towards a higher surface density ratio: at Eo = lOg cm"^ 
the relative change in outer semi-major axis, rj « 0.3, but for Eo = 
50g cm"^, this increases to = 0.79;^q24. Similarly, the mass of the 
outer-most body also increases significantly as the surface density 
and the surface density ratio are increased. 

In the extreme case that the number of planetesimals in the 
outer zone goes to zero, then the outward migration of the embryos 
from the inner zone is heavily suppressed, confirming the require- 



ment that there be an external population of some kind to enable 
scattering-driven migration to take place. 

We note at this point that as the surface density is changed 
between simulations, the initial number of embryos also varies (be- 
cause the isolation mass and semi-major axis separation are func- 
tions of surface density - see Eqn|2}, changing from 80 bodies when 
E = 5 to 25 bodies when E = 50. This change in the number of 
bodies does itself have consequences for the scattering rate (see 
g3.4.1l (. so we re-ran the simulations with the initial number of em- 
bryos artificially constrained to be a constant (56), and found that 
the change in the initial number of bodies (e.g. going from 25 to 56 
bodies at E = 50) caused only a small modification to the migra- 
tion and mass growth compared to that resulting from the overall 
surface density variation. 

Previous studies bv'Guede s et alj j200 8f) which looked at the 
evolution of embryos after the runaway-growth phase used a shal- 
lower surface density profile (E oc than that use d in the ma- 
jority of this work (E oc a 

-3/2). iThebault et al] ( |2009|) established 
that such a shallower profile did not help to extend the inner zone 
which is suitable for runaway growth. For the sake of completion, 
we investigate in set D whether a change in the surface density 
slope would have any effect on the outward migration in our n- 
body simulations. Changing the surface density profile obviously 
changes the relative masses at given points, so we set up simulation 
set D so that the absolute normalisations at lAC/ were different. 
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Figure 1. Top: Semi-major axis evolution in an example of the fiducial 
MMNM simulation. The embryos which initially have a < 0.7 AU are plot- 
ted in black, the planetesimals which are initially in the outer region are 
grey. For the outer-most embryo the masses (in units of M^) and relative 
change in semi-major axes (;; = (aouierjinal - aouter.imiial) I aouterjmtial) are 
labeled at ICf'yrs & lO'yrs. We see that significant embryo migration out 
to ~ 1 A [/ is observed. Bottom: The value of as a function of time, av- 
eraged over 10 simulations, plotting the median value, as well as the lower 
and upper quartile bounds. 



changing them so that the total mass in the inner zone for each 
simulation was identical (~ 2Mq). We found as expected that shal- 
lower slopes, which effectively increase the mass ratio between the 
outer and inner zones, increase the degree of outward migration, so 
that (at 10* vr.s), the relative outward migration, rj = 1.04+[j3^ for 
a completely flat (E oc a") simulation (compared to ;/ = 0.29 for 



3.3 Position of tlie Zone Boundary 

The simulations thus far have assumed that the boundary, a„-„, 
between the inner and outer zones lies at 0.7 AC, as appropriate 
for the nominal case of a Cen-B. However, the position of a„-;, 
is expected to vary s ignificantly for a number of reasons: (i) In 
iThebault et ^ ( l2009t) , Q.l AU was in fact the position that would 
be expected for a rather massive disk in which embryo formation 
extended to rather greater radii than would expected in a MMNM 
disk; (ii) Variations in inclination (see lXie & Zhoij2009l) have been 
found to change the size of the region over which runaway growth 
is suppressed; and (iii) Difl'ering relative binary parameters (mass 
ratio, separation and eccentricity) would all serve to change the lo- 
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Figure 2. Variation in scattering driven migration (top, = (aourer./ma/ ~ 
loiitir.imtial) I <^o\aer,inttial) and mass (bottom) due to changes in surface den- 
sity. The plots are (in order dark-to-light), ^outer I'^inner = 0; "^outer/dinner = 

1; "Zouier /dinner = 2; dgui^r /dinner = 3. Increases in surface density above 
the MMNM(10gcm-2) cause significantly increased migration. Plotted are 
the median values at 10^ yrs (10 simulations for each point), with the error 
bars giving the upper and lower quartiles. 



cation of a„.„. We therefore generalise our investigation to look at 
how our results will scale for other systems in which a„ii hes at a 
different distance from the central star. 

We note that as we keep the total mass and extent of our disk 
constant, moving a„„ in (out) has the effect of decreasing (increas- 
ing) the total number of isolation mass embryos that are simulated 
within the inner zone. We thus always calculate the 77-migration 
parameter by focusing on the outer-most embryo at any given time. 

In set E we look at the degree of outward migration that re- 
sults for disks in which Ocru varies between 0.4 and l.OAU. In ^3l\ 
we found that for our fiducial MMNM with a„i, = 0.7 At/, that the 
average relative change in semi-major axis of the outer-most em- 
bryo was T] = 0.3. Varying a„ii in Fig[3]reveals that at 10* years 
77 is largely unaffected by changes in a„ii' there is some scatter in 
the results obtained, but there is no obvious trend in the degree of 
migration, suggesting that the relative change in semi-major axis of 
the outermost embryo will effectively be independent of the precise 
location of a^.,-„. However, we note that the mass of the outer-most 
embryo declines as a function of as the embryos collide with 
a smaller number of planetesimals at larger semi-major axes. 

When we extend the simulations on to 10' years, as is depicted 
in the lower panel of Fig [3] we find that the situation is perhaps 
slightly more complicated than initially thought, as the value of rj 
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Figure 3. Variation in relative migration, r] = {aouitr.fiiml ~ 
aouterjmtial)laouter,inmal (soIid line, left-hand axis) and mass (dashed line, 
right-hand axis) of the outer-most embryo as a function of the boundary 
radius, acrit, between inner and outer zone. Each point gives the median 
value from 10 simulations, with the error bars giving the upper and lower 
quartiles. Top: lO^yrs - The relative migration is approximately constant at 
~ 0.3 suggesting that the relative change in semi-major axis of the outer- 
most embryo will be effectively independent of the precise initial value of 
acrit- The mass of this outer-most embryo shows a gentle decline as a func- 
tion of o„-„. Bottom: \(Pyrs - The variation in the relative migration, is 
large but again does not show an obvious trend as a function of Ucrit- 



at this time is now seen to vary more widely with a^rn . One could 
argue that there is perhaps an approximate peak in relative migra- 
tion for flcrir ~ 0.7 AU, but the variation in simulation results is so 
large (especially for low a„-;,) that it seems more sensible to take 
the view that the results are consistent with the relative migration 
being approximately constant at 77 ~ 0.6. 

Finally, we now consider what happens if we dispense with 
Ucrit completely, and simply consider a MMNM disk in which 
growth to isolation has not been suppressed for any a, i.e. the situ- 
ation as might be envisaged around a single (non-binary) star. We 
would expect it to have a chain of isolation mass embryos extending 
across the entire terrestrial region. We model such a disk, simulat- 
ing the evolution of a chain of embryos initially spread between 0.2 
and 1 .5 A (/. As might be anticipated, we find that there is a consid- 
erable degree of orbit crossing, scattering and collision amongst the 
embryos. 

If we consider just those embryos which start with a < OJAU 
we find that after lO^yrs (lO' yrs), the maximum semi-major axis 
to which this material has scattered / migrated, is such that rj = 
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Figure 4. Variation in the semi-major axis (black line, left-hand axis) and 
mass (grey line, right-hand axis) as a result of variations in the initial eccen- 
tricity of the embryo population. Results are measured att = 10^ yrs. It can 
be seen that as long as the eccentricity is "reasonable", i.e. less than 0.5, 
than there is no real trend for the maximum semi-major axis of the outer 
planet to vary with eccentricity. However, there does seem to be a trend for 
the mass to decrease with eccentricity, such that M x - log(e). 



0.29^15 13 (77 = 0.83^[5^;|). Thus there is a rather similar degree of 
"mixing" of material from a < acru to a > acru as was seen in set A 
in which the initial embryo growth beyond flrnv was suppressed. But 
it should be noted that these simulations with no a^.,.,, show much 
more variability than those in set A, with a third of the simulations 
showing no spread at all beyond 0.7 At/, where-as all of the sim- 
ulations in set A migrated to at least I.OAt/ by \QP yrs (with the 
average being 1.25 At/, 77 = 0.8). 



3.4 Varying the Composition of the Inner System 

3.4.1 Number and eccentricity of the inner embryo population. 

In set F we vary the initial content of the inner zone, keeping the 
total mass and surface density profile the same as in the fiducial 
case, but we now vary the number and individual masses of the 
embryos, and then examine the effect of varying their initial eccen- 
tricity. If the number of embryos in the inner zone is substantially 
lower than the 56 which results from assuming Aa = Iru (i.e. if 
we make K > 7), then this would be more consistent with a set of 
initial conditions occurring a substantial time period after the gas 
disk has dissipated and the subsequent stage of chaotic collision 
has proceeded for some time, resulting in more widely spaced (and 
more massive) embryos. 

As the initial number of embryos is decreased to 25, the 
amount of outward migration is decreased, with the outer-most 
embryo at 10'' yrs having migrated such that rj = 0.21^°Ug but the 
mass of the outer-most body increases to ~ 0.3Me, I.e. larger bod- 
ies have less chance of being scattered / migrating outwards. By 
the time the initial number of embryos is lowered to 10, the em- 
bryos are_so_widelyseparated, and the interaction timescales so 
long jZhou et al]|2007l) . that the bodies have little chance of being 
scattered away from their initial locations. 

In a scenario such as this (in which the embryos are more 
widely spaced as a result of an assumed period of scattering), we 
would expect that the eccentricities of the embryos would be sub- 
stantially higher, approa ching e ~ 0.2 a s the number of bodies 
drops into single figures ( IChambersll200lh . We investigate this, but 
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find from further simulation work (see Fig|4l( tliat tfie initial eccen- 
tricity is largely unimportant in determining the degree of migra- 
tion, as long as it is not above an improbably large value (~ 0.5). 
Interestingly however, we find that the mass of the outer object is 
influenced by the initial eccentricity, with initial eccentricities of 
e = O.OI (e = 0.2) resulting in masses for the migrating embryos of 
0.2Me (~ 0.1 M^), I.e. giving an approximate linear fit: 



Momr,Final{^=W yfS) 



0.08 - 0.06 log,o(e,„,M), 



(4) 



where we stress that this result is applicable to measurements taken 
after 10* years of the evolution of simulations of a IQgcnT^ sur- 
face density disk, in which isolation-mass embryos were initially 
confined to a region a < 0.7 At/. 



3.4.2 Introduction of Planetesimals to the Inner System. 

It is possible that at the time of disk dissipation, the inner zone may 
contain a mixed population of massive embryos and less-massive 
planetesimals that have not yet been accreted. In set G we vary the 
fraction of the inner zone which is composed of embryos, varying 
the fraction, /, from the the value of 1 used in previous sets down 
to 0.1. 

If we implement this variation by keeping the surface density, 
mass and number of embryos constant but add planetesimals to 
the inner zone (and also add planetesimals to the outer zone pro- 
portionately), then we have increased the overall surface density, 
and qualitatively the results are as might be expected from 33.21 
with both migration rates and growth rates increasing. We find at 
\(f yrs, that for / = 0.1 (i.e. meaning the total surface density is 
ten times higher), = 1.4;^g35. 

If we implement the variation by keeping the overall surface 
density constant, requiring the addition of planetesimal mass to the 
inner zone and the subtraction of embryonic mass (and leave the 
outer zone unaltered), then we find little difference between the 
/ = 0. 1 and / = 1 cases, with both having 77 ~ 0.3 at \& yrs. The 
planetesimals do not seem to have any significant damping eff^ect. 



3.5 Simulation Approximations 

3. 5. 1 Planetesimals 

Ideally we would have liked to perform our simulations using a 
vast number of very small planetesimals in the outer region, fully 
interacting with the embryo population as they orbited the central 
star under the influence of the perturbing binary. However, for the 
sake of computational efficiency and practicality we simulated only 
the planetary embryos as fully interacting particles, modelling the 
planetesimals as test particles with an associated mass. This means 
that gravitational embryo-embryo and embryo-planetesimal inter- 
actions took place in the simulations, but planetesimal-planetesimal 
interactions were neglected. We also used only a relatively small 
number of massive planetesimals (~ 100). In simulation sets H 
and I, we explicitly investigate the implications of these approxi- 
mations by conducting a limited number of simulations in which 
the number of planetesimals is increased and / or the planetesimals 
are modelled as fully-interacting n-body particles. 

Our fiducial set A simulations in which the mass of plan- 
etesimals in the outer disk is distributed amongst 100 planetesi- 
mals, found that the outermost embryo migrated outwards such that 
rj = 0.29 after IQf'yrs. When we distribute the planetesimal mass 
between 1,000 bodies (set H), we find that the outermost body has 
migrated by the same amount (again 77 = 0.31^1555*) in the same 
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Figure 5. We compare the evolution of the ;7-parameter (see Eqn[3) in sim- 
ulation set J (black line, includes binary perturber) with the results from set 
A (grey line, no binary perturber). The behaviour of the simulations is al- 
most identical, suggesting that the binary companion does not significantly 
change the outward migration. 



time period, suggesting that the degree and rate of migration are 
essentially unaffected when using planetesimal numbers ^ 100. 

Additionally, we find that it does not seem to make a signifi- 
cant difference to the migration history if we use a fully-interacting 
model for the planetesimals, rather than modelling them as mas- 
sive test particles. Switching between set A (no planetesimal- 
planetesimal scattering) and set I (planetesimal-planetesimal scat- 
tering allowed) resulted in the average maximum migration dis- 
tance of the outermost embryo being almost identical (0.89AI/). 
We therefore feel confident that our approximation of using a rel- 
atively small number of massive test-particles provides sufficient 
accuracy to allow us some insight into the behaviour of such sys- 
tems. 



3.5.2 Binary Companion 

As a final point, we turn to consider in greater detail whether it was 
reasonable to have ignored the presence of the binary perturber in 
our simulations. We now explicitly include the binary companion 
in our simulation set J, starting the simulations with a I.IMq mass 
companion at 18 At/ with e = 0.5 (i.e. using the parameters of the 
a Cen system, assuming the embryos are in orbit around a Cen- 
B). We also slightly alter the initial conditions for the embryo and 
planetesimal disk, taking into account the likely forced eccentricity 
due to the binary object. Using the standard first order expansion in 
the d isturbing function (Heppenheimer 1978; Murra v & Dermoid 
l2OO0l) . we find that e/orrerf ~ 0.035a. Other than this, the simula- 
tions are kept the same as Set A, with the same number of embryos 
and planetesimals, and the planetesimals performing the same scat- 
tering role. 

A comparison between the results of this simulation set (J) 
and the results of our fiducial set (A) is provided in Fig.|5] demon- 
strating that there is very little difference between the two results, 
suggesting that the binary companion is not important in dictating 
the degree of scattering-driven migration. 

Investigating further, we plot in Fig[6]the semi-major axis of 
the outermost embryo, demonstrating that its median value, as well 
as the upper and lower quartile spreads are essentially unchanged 
by the presence of the binary companion. The mass of the outer- 
most body in Fig[7]is also approximately the same for the more real- 
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Figure 6. Variation in average maximum embryo semi-major axis for sim- 
ulations with and without a binary perturber Green plots include the per- 
turber and a realistic forced eccentricity. Red plots do not have a perturber 
present in the simulation. The top panel displays the median semi-major 
axis, whilst the bottom panel also displays the evolution of the upper and 
lower quartile values from the 10 simulations conducted. The migration dis- 
tance is essentially unchanged by the presence of the perturber. 
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Figure 7. Variation in mass of the outer-most embryo for simulations with 
and without a binary perturber. Green plots include the perturber and a re- 
alistic forced eccentricity. Red plots do not have a perturber present in the 
simulation. The top panel displays the median mass, whilst the bottom panel 
also displays the evolution of the upper and lower quartile values from the 
10 simulations conducted. The median mass of the outer-most embryo is 
also unaltered by the presence of the binary, although it should be noted 
that the simulations without the binary have a larger spread in masses. 



istic simulation involving the binary companion plus initial forced 
eccentricity, but it should be noted that the simulations with the the 
binary companion (set J) have a much smaller spread in the upper 
and lower quartiles for the mass. 

Finally, we examine the eccentricities of the bodies in Fig [8] 
finding that both the eccentricity of the outer-most embryo and 
the average eccentricity of all the embryos evolve in rather sim- 
ilar manners: the simulations with a binary perturber start with a 
larger (forced) eccentricity, but over time this difference becomes 
of less significance as the embryos in all simulations start to stir 
themselves to similar levels. We note that in 33.4. 11 the accreted 
mass was seen to reduce slightly for simulations with higher eccen- 
tricity: if we substitute the comparatively small forced eccentricity 



Figure 8. We plot the evolution in the eccentricity of the embryos in simu- 
lation sets A and J, i.e. we compare the eccentricity evolution excluding and 
including the presence of the perturbing binary. In the top panel we show 
the average eccentricity evolution of the entire embryo population, whilst in 
the bottom panel we focus on the eccentricity of the outer-most embryo. In 
both panels the thin lines give the results of set A in which no binary com- 
panion is present , whilst the thick line gives the results for set J in which 
the binary companion is explicitly simulated. We find that the behaviour 
of both measures is rather similar, with the initial forced eccentricity being 
obvious in both panels, raising the eccentricities above the cases with no 
binary perturber. However, in both cases the difference between set A and 
set J diminishes with time as the embryos stir themselves and increase their 
eccentricity. 

(0.035) into Eqn|4] then we find that the expected reduction would 
be of only a few percent. 

We take this opportunity to reiterate that, even though the bi- 
nary companion was omitted from the majority of the n-body simu- 
lations, its suppression effect on runaway growth has of course been 
included via and the attendant lack of embryos for a > a„if 



4 CONCLUSION 

We have considered the formation and migration of protoplanetary 
embryos in disks around the stars in binary systems, in which the 
initial stages of embryo formation have been curtailed beyond some 
critical semi-major axis, a^rif We find that: 

• For a fiducial MMNM disk, over a period of lO^yrs (10'>Ti), 
systems will on average exhibit outward migration such that they 
contain embryos which have migrated from 0.7 At/ to 0.9 AU 
(1.2AU), i.e. they have experienced a relative change in semi-major 
axis of t] = 0.3 (0.8), with the average outer-most body having a 
mass of 0.2Me (0.4Ma5) (see Fig[T](. 

• The migration starts to plateau by the 10' yr mark in the sim- 
ulations, suggesting that little significant further migration is to be 
expected beyond this point (see Fig[TJ. 

• Higher surface density disks increase the rate and degree of 
migration as well as increasing the mass of the migrating body (see 
FigO. 

• The position of the zone boundary (a„-i,, outside which run- 
away growth is suppressed), is unimportant in deciding the degree 
of relative migration, zy, that takes place in the system: whilst there 
is significant scatter in the results obtained, they are consistent with 
T] being independent of a„-„ at a given time (see FigO. 

• The presence of the binary companion in the system, whilst 
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important in initially determining the cut-off radius for embryo 
growth, does not significantly impact the subsequent outward mi- 
gration rate of the embryos, nor their final masses (see Figs|5]|6]and 

m- 

We suggest that tight binary systems, where runaway embryo 
growth is initially confined to a < acru, may very generally expe- 
rience outward migration that is sufficient to move these planetary 
embryos to regions a > Ocrit- This allows the possibility that the 
habitable zone can be populated with detectable, terrestrial mass 
planets, even if a„» is significantly interior to the HZ. 
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